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导航 卫星 姿态 和 光 压 模型 研究 进展 
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摘要 : 
统 数 据 联合 处 理 时 ， 


黄 E, RANT, 


科学 院 上 海天 文 台 ， 


不 同 导 航 系统 中 卫星 的 形状 、 材 料 、 轨 道 
要 获得 


上 海 200030; 2. 中 


RAKHE, 


国 


设计 和 姿态 控 


高 精度 的 卫星 轨道 ， 需 要 区 别 对 待 不 


af 


By 1,2 


者 1 


科学 院 大 学 ， 北 京 100049) 


制 模式 等 都 不 相同 ， 在 进行 多 系 


司 类 型 卫星 的 姿态 模型 和 六 


GH 模 


型 。 作 为 导航 卫星 所 受 量 级 最 大 的 非 保守 力 ， 国 内 外 许多 研究 是 围绕 光 压 模型 的 建立 与 精 化 等 方 
面 来 开展 ， 以 便 建 立 最 合适 的 模型 ， 从 而 获取 更 高 的 定 轨 精 度 。 目 前 的 研究 表明 ， 新 的 Galileo 
以 及 北斗 卫星 采用 与 之 前 不 同 的 姿态 控制 模式 ， 适 用 的 光 压 模型 也 不 同 。 总 结 了 目前 比较 通用 的 
导航 卫星 姿态 模型 ， 尤 其 是 官方 最 新 公布 的 北斗 卫星 的 姿态 模型 ， 详 细 介绍 了 卫星 光 压 模型 研究 
进展 和 近 几 年 来 新 发 射 卫星 光 压 模型 以 及 其 先 验 模型 的 构建 方法 ， 并 分 析 了 各 个 模型 的 优 缺 点 
和 适应 性 。 通 过 对 目前 导航 卫星 姿态 以 及 光 压 模型 的 分 析 ， 提 出 在 更 高 精度 的 定 轨 任 务 中 有 待 改 
进 的 问题 ， 为 导航 其 是 新 发 射 导航 卫星 精密 定 轨 中 的 光 压 和 姿态 模型 的 选择 和 建立 提供 
参考 。 
X 键 i: 导航 卫星 ， 太 阳光 压 ， 卫 星 姿 态 ， 精 密 卫 星 定 轨 
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1 引 言 


欧盟 Galileo 和 我 国 北斗 导航 系统 (Beidou Navigation Satellite System, BDS) 作为 新 一 


代 的 全 球 导 航 卫 星系 统 (Global Navigation Satellite System, GNSS)， 目 前 已 经 提 


务 ， 起 步 更 早 的 GPS 和 GLONASS 也 在 i 


统 的 出 现 致力 于 提供 
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步 发 展 其 
全 球 覆 盖 、 全 天 候 和 高 精度 的 导 
实现 高 精度 导航 定位 服务 的 前 提 ， 精 密 卫星 定 轨 和 轨 
导航 卫星 受到 的 非 保守 力主 要 有 太 [ 
辐射 压 摄 动力 、 热 辐射 和 天 线 电磁 辐射 摄 动力 等 ， 
材料 属性 等 影响 "”。 


这 些 力主 要 


pay 


卫星 以 实现 现代 化 。 
航 定 位 授时 服务 。 精 密 的 卫 


供 全 球 服 


多 种 卫星 导航 系 


da 


预报 需要 
光 压 摄 动力 (solar radiation 


pals SRP), 
星 姿 态 、 几 何 结 


星 轨 道 是 


要 高 精度 动力 学 模型 的 支持 。 


地 球 反 照 
构 和 卫星 
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SRP 是 太阳 光子 撞击 物体 表面 产生 的 作用 力 ， 作 为 导航 卫星 所 受 量 级 最 大 的 非 保 守 力 ， 
SRP 误差 会 对 积分 得 到 的 卫星 速度 和 位 置 产生 累积 效应 ， 严 重 影响 精密 定 轨 的 精度 ”。 目 
前 ， 光 压 模 型 及 其 精 化 是 众多 国际 GNSS 服务 (the international GNSS service, IGS) 分 析 
中 心 的 热门 研究 内 容 ， 研 究 结果 表明 ， 不 同 的 SRP 建 模 方法 会 导致 不 同 IGS 分 析 中 心 的 轨 
道 产品 之 间 存 在 差别 "“。SRP 的 影响 因素 除了 总 太阳 辐 照度 (total solar irradiance, TSI), 
卫星 与 太阳 的 位 置 关 系 ， 还 有 卫星 本 身 的 姿态 、 面 质 比 、 形 状 、 材 料 等 。 姿态 控 制 的 精度 会 
直接 影响 到 光 压 模型 精度 ， 姿 态 控制 规律 是 建立 光 压 模型 的 基础 。 

对 于 导航 卫星 ， 目 前 可 用 的 光 压 模型 较 多 ， 部 分 卫星 应 该 采用 的 最 优 模型 以 及 相应 的 先 
验 模型 不 是 非常 明确 ; 对 于 同一 类 型 卫星 ， 在 轨 运 行 的 不 同时 期 (地 影 期 、 非 地 影 期 等 ) 其 
适用 的 光 压 模型 不 完全 一 致 ， 目 前 采用 异 构 星座 的 BDS 的 不 同类 型 卫星 光 压 模型 的 适用 性 
需要 进一步 验证 ， 最 新 发 射 的 GPS II、GLONASS-K 和 BDS-3 卫星 的 姿态 和 光 压 模型 也 
需要 验证 。 本 文 主要 介绍 了 各 导航 系统 的 姿态 模型 和 目前 常用 的 太阳 辐射 压 模 型 ， 分 析 了 各 
模型 的 适用 性 ， 为 导航 卫星 的 太阳 光 压 建 模 及 模型 优化 提供 借鉴 。 


Ch 


2 导航 卫星 姿态 控制 


2.1 卫星 姿态 控制 规律 
3J 星 姿态 是 指 卫星 各 轴 在 轨道 上 所 处 的 空间 指向 状态 。 导 航 卫 星 为 实现 正常 运行 ， 其 天 
线 要 指向 地 心 ， 太 阳 帆 板 平 面 要 垂直 于 卫星 到 太阳 的 连 线 。 在 名 义 姿态 下 ， 卫 星 姿态 满 足以 
上 两 点 。 为 保持 名 义 姿态 ， 卫 星 姿态 控制 系统 需要 不 断 进行 调整 。 
卫星 姿态 对 精密 定 轨 的 影响 体现 在 两 个 方面 : (1) 卫星 天 线 的 相位 中 心 偏差 和 变化 等 几 
何 误差 的 改正 ; (2) 卫星 姿态 引起 受 照 面 的 变化 ， 进 而 导致 太阳 光 压 摄 动力 、 地 球 反照 辐射 
压 摄 动力 等 发 生变 化 。 因 此 ， 为 了 构建 高 精度 光 压 模型 并 实现 导航 卫星 精密 定 轨 ， 首 先 要 获 
得 准确 的 卫星 姿态 呈 。 精 密 定 轨 时 卫星 姿态 主要 用 来 描述 星 固 坐 标 系 与 惯性 坐标 系 之 间 的 转 
换 关 系 。 
GNSS 卫星 星 固 系 的 原点 置 于 质心 ，Y 轴 为 太阳 帆 板 的 旋转 轴 方向 ，2 轴 为 治 卫星 工 
波段 天 线 发 射 信号 方向 ，X 轴 与 以 上 两 轴 正 交 构 成 右手 坐标 系 。 图 1 给 出 星 固 系 及 卫星 、 
地 球 、 太 阳 三 者 相互 的 位 置 关系 ”。 其 中 ，6 为 太阳 高 度 角 ， 表 示 太 阳 和 地 心 连 线 与 轨道 面 
之 间 的 夹 角 ; «© 为 从 卫星 上 看 地 心 与 日 心 之 间 的 角 距 ; / 为 轨道 角 ， 表 示 在 轨道 面 内 远 日 点 
和 地 心 连 线 与 卫星 和 地 心 连 线 之 间 的 夹 角 。X 轴 反 映 卫 星 深 动 方向 ，Y 轴 反 映 卫 星 俯仰 方 
H, Z 轴 反 映 卫星 偏 航 方 向 。 导 航 卫 星 三 轴 姿 态 稳定 模块 可 以 控制 星 固 系 2 轴 指 向 误差 不 
大 于 0.1*， 滚 动 和 俯仰 的 影响 可 以 忽略 ”。 导 航 卫 星 的 姿态 控制 主要 是 偏 航 控制 。 图 2 和 图 
3 为 中 国 空间 技术 研究 院 (China Academy of Space Technology, CAST) 和 中 国 科学 院 上 海 
微小 卫星 创新 研究 院 (Shanghai Engineering Center for Microsatellites, SECM) 研制 的 中 圆 
轨道 (medium Earth orbit, MEO) 卫星 示意 图 。 


A 
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1 卫星、 地球、 太阳 位 置 关 系 及 星 固 坐标 系 


2 ”中 国 空间 技术 研究 院 MEO 卫星 示意 图 3 ”上 海 微小 卫星 工程 中 心 研制 MEO 卫星 示意 图 


名 义 偏 航 角 为 星 固 系 X 轴 正 向 和 卫星 速度 方向 形成 的 夹 角 ， 按 下 式 计算 : 


W = ATAN2(—tan@,sinp) , (1) 
W = ATAN2(tan $, — sin y) , (2) 


其 中 ，ATAN2 等 于 arctan(a/b) 的 计算 值 ， 返 回 值 范围 为 [一 180?, +180°], 2, u 与 图 1 中 表 
示 一 致 。 对 于 不 同 卫 星 ， 当 轴 X 正 向 指向 太阳 时 ， 按 式 (1) 计算 ， 采 用 这 种 模式 的 卫星 有 
GPS Block IIF、GLONASS-M、 北 斗 倾斜 同步 轨道 (inclined geosynchronous satellite orbit, 
IGSO) 和 MEO 等 ; 当 轴 X 正 向 背 向 太阳 时 ， 按 式 (2) 计算 ， 采 用 这 种 模式 的 卫星 有 GPS 
Block IIR. Galileo 等 5"。 两 种 模式 差别 体现 在 名 义 轴 X 与 Y 不 同 。 由 于 姿态 控制 系统 动 
量 轮 调整 速率 有 限 ， 目 前 GNSS 卫星 姿态 控制 策略 主要 有 两 种 ， 一 种 是 在 太阳 高 度 角 小 于 
某 个 阔 值 时 进行 轨道 机 动 ， 另 一 种 是 在 动量 轮 调整 速率 大 于 一 定 值 时 开始 轨道 机 动 ”。 


334 


天 文学 进展 


将 式 (1) 与 (2) 对 时 间 求 偏 导 ， 可 以 得 到 偏 航 角 速率 : 


KF, 为 卫星 平均 运动 角 速 率 。 在 u NO. 180°R, U BKA | 


点 时 ， 卫 星 偏 航速 率 有 可 


2.2 ” 现 有 导航 卫星 姿态 模型 
目前 各 GNSS 系统 主要 采用 3 种 姿态 控制 模式 : 名 义 偏 航 (Smt h 


2s L tan 8 cos u 


pi 


sin? u + tan? 8 


能 超过 动量 轮 调整 速率 限 值 ， 产生 轨道 机 动 
要 足够 的 光照 ， 当 太阳 高 度 角 较 小 时 ， 某 些 卫 星 还 会 产生 地 影 机 动 。 


由 于 


41 卷 


(3) 


Ey 即 在 近日 点 、 运 日 
卫星 光敏 感 器 需 


D) 模式 、 机 动 偏 航 


控制 (连续 控制 偏 航 ) 模式 和 零 偏 控制 模式 。 当 偏 航 角 随 太 阳 、 地 球 和 卫星 的 位 置 变化 而 变 
化 ， 使 卫星 保持 名 义 姿态 时 称 为 名 义 偏 航模 式 ， 为 了 保证 太阳 帆 板 法 向 与 太阳 光线 之 间 夹 角 


不 大 于 设 定 值 ， 侦 航 角 固 


近 以 设计 的 姿态 变化 或 以 最 大 角速度 变化 以 实现 目标 偏 航 角 ， 称 为 机 动 


定 为 0"， 称 为 零 偏 控制 模式 ， 对 于 地 影 期 卫星 ， 在 午夜 、 正 午 附 


ee 
此 类 卫星 地 影 期 的 范围 大 ， 所 以 不 存在 
第 一 种 ，BLOCKIIF、BLOCKI VIA 卫星 的 轨道 机 动 模式 为 第 二 种 。 
动 与 午夜 机 动 的 模型 基本 一 致 ， 当 太阳 高 度 角 小 于 设计 值 # 
轮 按 照 最 大 角速度 旋转 ， 直 到 实际 姿态 与 名 义 姿态 一 致 。B 


根据 出 地 影 时 的 偏 航 角 计 算 回 
( 5 min LAA), 可 以 忽略 地 影 
前 在 轨 运 行 的 GLONASS 卫星 绝 大 多 数 为 GCLONASS-M， 后 续 
在 设计 上 与 GLONASS-M 主要 不 同 的 是 采 /) 


K1 和 GLONASS-K2 TÆ, 


星 地 影 机 动 有 所 差别 ， 对 了 
由 于 出 地 影 时 的 卫星 姿态 与 名 义 姿态 并 不 一 致 ， 因 此 存在 地 影 恢复 期 ， 对 于 BLOCKIIF, 
定 角 速度 进行 地 影 机 动 ，BLOCKIIF 卫星 的 地 影 恢 复 期 较 短 


偏 控制 模式 相 比 ， 机 动 偏 航 控制 模式 显著 解决 了 服务 性 能 下 降 的 问题 


E 午 机 动 和 午夜 机 动 ， 男 一 种 有 
午夜 机 动 。 其 中 ，BlockILR 卫星 的 轨道 机 动 模式 为 
各 类 卫星 的 正午 机 
达到 开始 机 动 的 轨道 角 时 ， 动 量 
LOCK IIF #1 BLOCK II/ILA E 


局 航 控 制 模 式 。 与 零 


E 午 机 动 、 地 影 机 动 ， 


F BLOCKIL/IA 卫星 ， 动 量 轮 以 最 大 角速度 旋转 ， 直 到 出 地 影 ， 


恢复 期 的 影响 “” 。 


发 射 的 GLONASS- 


窝 板 的 非 密闭 星体 设计 ” 


> GLONASS-M 卫星 轨道 机 动 模式 为 地 影 书 


|B] < 14.2° 时 ， 进 入 地 影 机 动 ， 卫 星 将 以 约 0.25 (*):s-! 的 速率 调整 ， 


出 地 影 时 的 名 义 姿态 ， 之 后 


Ñ 
制 规律 ， 使 得 卫星 姿态 以 了 


< 


Ay 


T | 


卫星 偏 航 角 不 变 ， 直 到 出 地 影 ， 在 正午 机 动 时 ， 
某 一 时 刻 提 育 前 以 最 大 调整 速率 进行 姿态 控 制 ， 一 直 保 持 到 经 过 正午 点 后 卫星 满足 姿态 控 


几 动 和 正午 机 动 。 当 


j 了 一 种 基于 蜂 


以 最 | 


决 的 速度 调整 到 


J 星 是 在 正午 点 


F 为 中 心 对 称 分布 ， 实 现 整个 轨道 机 动 对 月 


昌 户 影响 最 小 ”"。 目 前 


的 研究 表明 ，GLONASS-K 在 地 影 期 采用 与 GLONASS-M 卫星 相似 的 姿态 控 


对 于 Galileo, 2017 年 10 月 ， 为 了 更 好 地 完成 Galileo FAC AY AHH Mb 


A, 欧洲 全 球 卫星 


Hil Mes! 


导航 系统 管理 局 公开 了 包含 偏 航 姿态 模型 在 内 的 各 卫星 的 相关 参数 。 Galileo 卫星 在 发 生 轨 


道 机 动 时 采用 机 动 偏 航 控制 的 姿态 ， 


Y = ATAN2(Shy, Se), 


Sz, Sy, S 为 太阳 单位 向 量 


E, 


可 根据 太阳 高 度 角 8 和 轨道 角 A (以 近 


点 为 起 点 


IOV 卫星 在 轨道 机 动 期 间 的 姿态 模型 为 : 


(4) 
) 计算 得 到 ， 
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表达 式 为 : 
[S. Sy S.]" = (— sin u cos 8 — sin 8 — cosucosB)' ， 


Si,, 为 平滑 后 的 S,. FOC 卫星 在 轨道 机 动 期 间 的 姿态 模型 为 ; 


Y = 90° - SIGN(1, (p15) + [P (us) — 90° - SIGN (1, Y (us))] : cos (Z £e) ， (5) 


HEP, ps 为 开始 轨道 机 动 时 的 轨道 角 ， 亚 (1) 为 开始 轨道 机 动 时 的 名 义 姿态 角 ，SIGN (a,b) 
的 返回 值 取 5 符号 和 a 的 值 ， 如 ax 表示 轨道 机 动 的 最 长 时 间 ，k 为 轨道 角 的 平均 变化 速率 。 
Galileo 的 姿态 模型 保证 卫星 偏 航 角 速率 在 0.203 (%).s-! 以 内 ， 姿 态 机 动 起 止 的 位 置 关 于 午 
夜 (或 正午 ) 对 称 ， 保 证 在 正午 (或 午夜 ) 时 其 实际 偏 航 角 与 名 义 偏 航 角 一 致 (90°)”。 

对 于 BDS-2， 郭 靖 " 建立 了 IGSO 和 MEO 卫星 的 姿态 模型 : |8| < 4° 为 零 偏 模式 ， 
18| > 4° 为 名 义 偏 航模 式 。 在 此 基础 上 。Dai 等 人 ”” 建立 了 更 加 精确 的 姿态 模型 ， 名 义 偏 航 
转 零 偏 的 条 件 为 : |6| < 4° A |W] < 5° MH |B) << 4° H5 < < 20°, WW>0; 零 偏转 名 
义 偏 航 的 切换 条 件 为 : |6| > 4° 且 人 要 | 乏 5 BR |B] > 4° 且 5 < |W] < 20°, WW>0. 2016 
年 3 月 发 射 的 BDS-2 卫星 IGSO-6 (C 016, C13) 的 姿态 控制 模式 并 不 是 名 义 偏 航 和 零 偏 ， 
而 是 名 义 偏 航 和 机 动 偏 航 控制 ， 这 种 新 的 模式 与 伽利略 FOC 卫星 的 姿态 模型 非常 类 似 。 此 
外 ，MEO-6 (C015, C14) 在 2016 年 10 月 后 遵循 机 动 偏 航 控制 模式 ，IGSO-1 (C05, C06) 在 
2017 年 3 月 后 也 遵循 机 动 偏 航 控制 模式 ， 说 明 在 轨 卫 星 的 姿态 控制 模式 也 可 以 设置 。 目 前 
所 有 BDS 的 倾斜 同步 轨道 卫星 (geostationary Earth orbit, GEO) 卫星 均 采 用 零 偏 姿态 控制 
模式 ”。 对 于 CAST 的 BDS-3 卫星 ， 当 太阳 高 度 角 在 3° 之 间 ， 且 太阳 -地 球 矢 量 在 轨道 
面 内 投影 与 地 球 -卫星 连 线 之 间 的 夹 角 < 30° 时 ， 偏 航 角 近似 为 : 


g Í ATAN2(tanb, Sos) 8>0 
—ATAN2(tan ð, Sos) B<0 


， (6) 


` 
E 


Om = 3/57.3(rad) 》 
8 = bm (cos(0.596s))> 4 


ERP, Ss 为 轨道 系 太 阳 矢 量 的 和 Was, So. 为 轨道 系 太 阳 矢 量 的 2 方向 分 量 。 
SECM 研制 的 BDS-3 卫星 采用 的 是 如 下 偏 航 姿态 模型 ”: 


| ATAN2(—tanp,sinu) |B] > Bo 
y= 


ATAN2(— tan bo,sin u) 0< 8< bo ‘ (7) 
ATAN2(tan o, sin p) -B <B <0 


tH, 6 = 3。 是 设 定 阔 值 。 目 前 的 研究 表明 : 卫星 的 实际 姿态 与 目前 公布 的 姿态 模型 并 不 
完全 一 致 ， 在 太阳 高 度 角 较 低 时 ，BDS-3 卫星 的 姿态 变化 需要 进一步 研究 。 表 1 为 目前 常 
导航 系统 卫星 的 姿态 控制 模式 
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表 1 导航 卫星 姿态 控制 模式 


太 
SNL AE 适用 时 姿态 模型 
航 系 统 控制 模式 适用 时 段 姿态 模型 
BDS-2(GEO) Shi 所 有 时 段 偏 航 角 为 0 
yg (BL |Ø] > 5° R |B| > 4, esr ; 
BDS-2(1) CONN 5 |W] < 20°, WH > 0 时 EE) Ao 
FE Ape = |B) < 4°, |Y SOAS: 偏 航 角 为 0 
5° <|V| < 20°, YY > 0 H 
imme H223 , |u — 180°] > 5° 或 ak n 
BDS-3(CAST) 名 义 偏 航 iu] > 5° 时 MATHAIR (1) 计算 
DERE ses 当 |B| < 3°, |u — 180°] < 5° 时 RSIR (4) 计算 
pri 4 |B) < 3°, |p| < 5° 时 MIZ (4) 计算 
名 义 偏 航 = |8| > 3° 时 扁 航 角 按 式 (1) 计算 
BDS-3(SECM) sn anim 4 JB] < 3° 时 采用 |B] = 3° 的 理论 偏 航 姿态 
18| > 正午 或 子夜 机 动 的 限 值 ed eae al : 
GPS Block aii = gi EEN 
ae SOUR ”或 未 达到 开始 机 动 的 轨道 角 时 We i 
47 AE fe 速度 
机 动 偏 航 。 |6| < 正午 或 子夜 机 动 的 限 值 ， BAe ernie 
影 ~ fs [ 动 4 L 道 . HE 
(地 影 ) 达到 开始 机 动 的 轨道 角 时 影 机 动 存在 地 影 恢复 其 
i [3] > 24° 或 未 达到 开始 Fad 要 
Do lode i 名 义 偏 航 机 动 的 轨道 角 时 高 航 角 按 式 (2) 计算 
By blo |B| < 2.4*， 达 到 开始 按照 最 大 角速度 旋转 ， 
TaN 机 动 的 轨道 角 时 直到 达到 名 义 偏 航 角 
|B] > 4° 或 未 达到 开始 
正午 机 动 的 轨道 角 时 ; 
i o 偏 航 角 按 式 (1) 计算 
GPS Block IIF IB8| > 8” 或 未 达到 开始 
午夜 机 动 的 轨道 角 时 
榨 昭 最 y 度 
机 动 偏 航 |B] < 各 ， 达 到 开始 2 o ee 
(正午 ) 正午 机 动 的 轨道 角 时 直到 达到 名 义 偏 航 5 
aT 也 影 机 动 存在 地 影 恢 复 期 
机 动 偏 航 18| < 8"， 达 到 开始 根据 出 地 影 时 的 偏 航 角 
(地 影 ) 午夜 机 动 的 轨道 角 时 计算 固定 角速度 
|B] > 2° 或 未 达到 开始 
午 机 动 的 轨道 角 时 ， 
eg SAE 偏 航 角 按 式 (1) 计 算 
GLONASS |6| > 14.2° 或 未 达到 开始 
午夜 机 动 的 轨道 角 时 
机 动 偏 航 18| < 2， 达 到 开始 
(正午 ) 机 动 的 轨道 角 时 保证 偏 航 姿态 以 正午 点 
为 中 心 对 称 分 布 
SAY ee dh y PEE 
机 动 偏 航 18| < 4.2*， 达 到 me o 
Ae. 开始 机 动 的 轨道 VIA EJ D Fh 
(午夜 ) 于 始 机 动 的 轨道 角 时 的 名 义 偏 航 角 
ed ea 名 义 偏 航 |6| > 2° 3È |u — 180°| > 15° 时 偏 航 角 按 式 (1) 计算 
机 动 偏 航 |6| < 2°, |u — 180°| < 15° 时 偏 航 角 按 式 (4) 计算 
Grane 名 义 偏 航 ”|8| > 4.1° EÈ |u — 180°| > 10° 时 偏 航 角 按 式 (2) 计算 
机 动 偏 航 |6| < 4.1°, |u — 180°] < 10° 时 偏 航 角 按 式 (5) 计算 
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3 光 压 模型 


导航 卫星 姿态 模型 影响 导航 卫星 包括 太阳 光 压 在 内 的 非 保 守 力 ， 是 导航 卫星 精密 定 轨 
和 光 压 模型 建立 的 前 提 。 除 卫星 姿态 和 卫星 形状 、 材 料 等 本 身 的 物理 特性 外 ，SRP 还 会 受 
TSI 和 卫星 受 太阳 照射 部 分 的 面积 影响 ， 而 且 这 些 因 素 始终 在 变化 ， 都 应 该 作为 光 压 建 模 
的 考虑 因素 之 一 。 研 究 表明 ， 真 实 太阳 活动 的 周期 性 变化 幅 值 小 于 0.1%， 基 于 变化 的 TSI 
模型 的 GPS 定 轨 精 度 略 有 提高 ， 提 高 幅度 在 0.1.0.5 mm”. XF BDS 的 倾斜 同步 轨道 
IGSO 卫星 和 中 国 轨道 MEO 卫星 ， 提 高 幅度 在 0.5 mm 以 内 ; 对 于 BDS 的 地 球 静 止 轨道 
GEO 卫星 ， 提 高 幅度 可 达到 3 mm ”。 此 外 ， 由 于 导航 卫星 在 轨 和 运行 时 ， 卫 星 的 一 部 分 会 
被 地 球 (或 月 球 ) 遮挡 ， 地 影 因 子 ( 蚀 因 子 ) 用 来 表示 卫星 在 被 遮挡 时 太阳 可 见 部 分 的 视 面积 
相对 于 完整 太阳 的 视 面 积 的 比例 ， 因 此 还 需要 计算 卫星 的 地 影 因 子 以 精确 获得 SRP。 导 航 
卫星 SRP 模型 主要 可 以 分 为 三 种 : 分 析 型 、 经 验 型 模型 和 半分 析 型 ， 下 面 主 要 介绍 目前 常 
用 的 几 种 光 压 模型 。 
3.1 “分 析 型 光 压 模 型 
3.1.1 ROCK 系列 模型 
ROCK 模型 是 GPS 卫星 的 制造 商 Rockwell 公司 的 Porter 等 人 (1976) 根据 卫星 的 外 部 
攻 状 和 材料 等 设计 的 太阳 光 压 模型 。 该 模型 在 星 固 系 下 描述 ， 通 过 简化 卫星 的 结构 ， 根 据 
卫星 的 光 压 参数 、 地 影 因子 、 各 部 分 之 间 的 遮挡 关系 等 ， 建 立 了 光 压 模型 ”。 鉴 于 Block I 
和 Block I 卫星 在 形状 、 材 料 光学 特性 等 方面 的 差异 ， 分 别 建立 了 ROCK4 模型 (1985) 和 
ROCK42 模型 (1989)”” 
对 于 任意 平面 ， 该 模型 包括 二 个 力 分 量 ， 分 别 由 入 射 光 、 镜 面 反射 光 、 漫 反射 光 产 生 。 
其 中 法 向 分 量 垂直 于 表面 ， 它 由 入 射 光 的 法 向 分 量 加 上 镜面 反射 光 的 反 冲 产生 ， 可 表示 为 : 
机 一 一 (4B/c(L+Hz)cos20 , (8) 
其 中 , v IRIX, 为 镜面 反射 率 (v 和 的 范围 都 为 0~1)，4 为 平面 的 面积 ，0 为 入 射 
太阳 光线 与 平面 法 线形 成 的 角度 。 光 压 常数 是 太阳 功率 常数 巨 ( 约 1368 wm?) 除 以 光速 。 
切 向 分 量 与 平面 相 切 ， 远 离 太 阳 。 它 是 由 入 射 光 的 切 向 分 量 减 去 镜面 反射 光束 带 走 的 动量 而 
产生 的 ， 可 表示 为 : 


He 


Fıs = —(AE/c)(1 — pv) sin ð cos ， (9) 
漫 反 射 分 量 与 法 向 分 量 位 于 同一 方向 ， 可 表示 为 : 
Fip =- (AE/c)v(1 — p)cos6 . (10) 


对 于 圆柱 体 来 说 ，3 个 分 量 分 别 为 : 
Fn = —(AE/c)(1 + pv) cos? 6 
Fəs = —(AE/c)(1 — pv) sin 8 cos 0 (11) 


Fyn = —7(AB/o)u(1 — p) cos 8 
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其 中 ，4 为 过 圆柱 体 轴 的 截面 面积 。 

民 据 各 卫星 的 部 件 结构 和 光学 属性 参数 ， 将 所 有 辅助 角 用 太阳 和 +2 方向 的 夹 角 BO 
示 ， 从 而 实现 仅仅 以 B 为 自 变量 分 别 计算 各 部 分 在 星 固 坐标 系 轴 X 和 轴 2 方向 的 SRP 
DE RYH ROCK 模型 称 为 标准 ROCK 模型 (ROCK S 模型 ，Block I 对 应 模型 是 S10, 
Block II/IIA 是 S20)， 其 中 S10 模型 为 (单位 为 10-5 N): 


X = —4.34sin(B) + 0.10sin(2B + 1.1) — 0.05 cos(4B) + 0.06 


12 
Z = —4.34cos(B) + 0.17 sin(2B — 0.4) — 0.05 sin(4B) — 0.06 We 
S20 模 型 为 : 
X = —8.10sin(B) + 0.05 cos(2B) — 0.056 sin(4B + 1.4) + 0.07 as 
Z = —7.80 cos( B) + 0.024sin(2B — 0.8) — 0.047 sin(4B + 0.9) — 0.02 


此 后 ， 考 虑 到 Block II 卫星 的 升级 以 及 卫星 自身 的 热 辐 射 ， 建 立 了 新 的 热 ROCK 模型 
(ROCK T 模型 ，Block I 是 T10，Block II/IIA Æ T20, Block IIR 是 T30)”“。 升 级 后 的 
T10 模型 为 : 


X = —4.55sin(B) + 0.08 sin(2B + 0.9) — 0.06 cos(4B + 0.08) + 0.08 


14 
Z = —4.54sin(B) + 0.20 sin(2B — 0.43) — 0.03 cos(4B) ue 
T20 模型 为 : 
X = —8.96sin(B) + 0.16 sin(3B) + 0.10 cos(5B) — 0.07 sin(7B) (15) 
Z = —8.43 cos(B) 
T30 模型 为 : 
X = —11.000sin(B) — 0.20sin(3B) + 0.20 sin(5B) T 


Z = —11.30 cos( B) + 0.10 cos(3B) + 0.20 cos(5B) 


3.1.2 Box-Wing 模型 

鉴于 卫星 形状 较 复杂 ，Marshall 和 Luthcke™” (1992) FET —PAR (Box-Wing) 模型 
IF Topex/Posiden 卫星 的 精密 定 轨 。 在 给 定 各 表面 的 几何 、 物 理 参数 及 卫星 姿态 信息 的 
前 提 下 ， 可 以 计算 出 各 个 表面 产生 的 太阳 光 压 摄 动力 ， 将 各 个 面 的 摄 动力 登 加 就 得 到 总 的 光 
压 摄 动力 “。 该 模型 能 够 反映 卫星 受到 的 太阳 光 压 摄 动力 ， 属 于 分 析 模型 。 模 型 将 卫星 的 
部 件 简化 ， 将 卫星 星体 的 六 个 面 和 太阳 能 帆 板 的 两 个 面 分 别处 理 ， 作 用 在 其 中 一 个 平面 的 太 
阳光 压 可 写 为 “”: 


23] 


0 
amca e ja p)en + 2 (3 | poso ) ex , (17) 


其 中 ，p 为 反射 光子 分 数 (镜面 反射 系数 )，6 为 散射 光子 分 数 ( 漫 反射 系数 )， 另 外 还 有 一 
个 吸收 光子 分 数 a (吸收 系数 )， 而 且 a 十 p 十 6 = 1， 可 以 由 式 (8) 中 的 反射 率 和 镜面 反射 
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率 转换 得 到 ，a = 1 一 v, p= w, ô= vl- u); A 是 对 应 平面 的 表面 积 ，M 为 总 质量 ; 
So = 1367 W-m-?, JQARAPAAY TSI; c 为 光速 ，ep 和 en 为 方向 向 量 ， 分 别 指向 太阳 方向 
和 平面 法 向 ; 0 是 ep 与 en 之 间 的 夹 角 。 卫 星 太阳 能 帆 板 所 受 的 光 压 力 可 以 按 式 (17) 计算 ， 
对 于 卫星 本 体 ， 表 面 有 绝缘 材料 ， 根 据 兰 伯 特 定律 ， 卫 星 表面 吸收 的 能 量 会 以 热 的 形式 被 反 
射 回 去 ， 在 式 (17) 中 减 去 反射 的 光 压 力 ， 卫 星 的 本 体 所 受 的 光 压 力 可 以 按 下 式 计算 : 


AS, [ 2 
a 一 Say cos 0 Ke +ô) (ev F Sex) + 2p cos dex : (18) 


分 析 模 型 有 清晰 的 物理 意义 ， 便 于 理解 ， 但 模型 公式 一 般 相 对 比较 复杂 ， 而 且 这 类 模型 有 明 
显 的 缺点 : (1) 对 于 有 复杂 结构 及 产生 遮挡 的 卫星 来 说 ， 难 以 准确 建 模 ， 分 析 型 光 压 模型 仅 
能 描述 光 压 的 主要 部 分 ， 并 不 能 准确 计算 卫星 受到 的 光 压 力 ; (2) 卫星 在 轨 上 服役 时 间 较 长 ， 
一 般 有 十 几 年 ， 由 于 卫星 表面 材料 老化 ， 卫 星 在 轨 运 行 的 不 同 阶段 受到 的 光 压 力也 在 不 断 变 
W” 
3.2 ”经 验 型 光 压 模型 
3.2.1 Colombo 模型 

在 经 验 型 模型 方面 ，Colombo (1989) 发 现 GPS 卫星 的 广播 星 历 轨道 误差 存在 共振 现 
象 ， 可 以 采用 一 个 经 验 模型 来 吸收 该 误差 ， 从 而 较 好 地 消除 包含 地 球 重力 场 、 太 阳光 压 在 内 
的 各 种 摄 动力 误差 对 定 轨 的 影响 ""。 该 模型 可 以 写 为 两 个 共振 频率 为 0 和 no 的 傅 里 叶 分 量 
ZA: 


Aa,, = Az, Cos not + B,, Sin not + Cy, , (19) 


其 中 ,i = 1, 2, 3 时 分 别 对 应 三 个 方向 的 轨道 误差 ， 待 估 参 数 为 A,,、B,,、0C,,。 但 是 该 模 
型 在 用 于 改进 GPS 卫星 轨道 时 认为 地 球 重 力 场 是 主要 的 误差 源 ， 对 于 卫星 与 太阳 位 置 关 系 
并 没有 明确 的 相关 性 ， 因 此 该 模型 对 于 SRP 的 改进 并 不 完善 。 
3.2.2 ECOM 系列 模型 
3.2.2.1 原始 ECOM 模型 

瑞士 伯尔尼 大 学 的 Beutler 等 人 在 1994 年 根据 经 验 的 Colombo 模型 ， 采 用 欧洲 定 轨 
中 心 (the Center for Orbit Determination in Europe, CODE) 的 数据 开发 了 CODE 经 验 轨 
道 模型 (empirical CODE orbit model, ECOM)， 模 型 采用 的 DY B 坐标 系 将 光 压 分 解 为 三 
个 很 适合 光 压 建 模 的 方向 D, Y, Bo HAI 可 知 ，D 方向 为 卫星 到 太阳 方向 ，Y 方向 为 大 
阳 帆 板 旋转 轴 ，B 方向 与 以 上 两 个 方向 正 交 构成 右手 系 。ECOM 将 待 估 参 数 扩 展 到 9 个 
(DO. YO. BO 和 6 个 周期 项 参数 )， 并 结合 ROCK 模型 提供 先 验 光 压 信息 ， 获 得 了 较 优 的 
定 轨 精 度 ”。 具 体 模型 为 : 


QsRP = aRock +4ap:ep +tay:ey +4B'EB , (20) 


右边 第 一 项 为 ROCK HAHAE, ap ay. agp 为 D、Y、B 三 个 方向 的 SRP， 分 
别 为 : 
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ap = Do + De cosu + D, sinu 
ay = Yọ + Y. cosu + Y, sinu y (21) 
ag = Bo + Be cosu + B, sin u 
KP, u 为 升 交 点 角 距 ， 范 围 为 [0, 360] (注意 与 轨道 角 / 进行 区 分 )。 尽 管 9 参数 ECOM 
在 研制 之 时 尚未 发 射 Block IR 卫星 ， 但 随后 的 应 用 中 发 现 该 模型 同样 适用 于 Block IR E 
星 。 该 模型 自 1996 年 应 用 于 CODE 分 析 中 心 。9 参数 会 增加 计算 的 自由 度 ， 对 于 三 天 解 


来 说 ， 不 需要 解 算 全 部 9 个 参数 ， 而 且 轨 道 参数 之 间 以 及 与 


A 


HIR H 


明 ， 忽 略 ECOM 在 B 方向 的 参数 ， 而 加 入 5 个 随机 力 参 数 ， 能 显著 提高 定 轨 精度 ( 


5~ 15cm) 


3.2.2.2 ECOMI 模型 


Springer 等 人 (1999) ZIL, H 
此 提出 简化 的 ECOM 模型 


入 1 


iz, 因 


数 。ECOM1 对 于 GPS 卫星 可 获得 更 


于 参数 之 间 的 相关 性 


~ 


他 参数 (如 UT1-UTC) 之 间 


E， 当 弧 段 长 度 增 加 到 5 d 时 ， 相 关 性 降低 到 可 接受 的 水 平 。 进 一 步 的 下 


究 表 


E， 太 多 光 压 参数 会 在 其 他 估 值 9 


GA 


el 


研究 表明 ， 与 其 他 赂 


有 相关 性 ，ECOMI1 模型 避免 了 这 种 问题 
共有 5 个 参数 ， 具 体 模型 为 : 
QD = Do 
ay = Yo 
ap = Bo + 


的 ECOM 相 比 ， 采 


Block WA 卫星 开发 了 


参数 与 8 的 关系 ， 采 用 5.5 a 的 轨 


新 的 先 验 模型 


ples 


模型 提高 了 一 个 量 级 
2013 年 一 直 采 用 该 先 
行 精密 定 轨 ， 也 能 获 
ECOM1”” 


o 


研究 表明 ，ECOM1 适 月 
不 过 对 Galileo 和 QZSS 1 


误差 。 这 种 现象 主要 原因 是 : Galileo 和 QZSS 卫星 主体 是 长 方 体 ， 卫 星 收 到 的 太阳 光 压 摄 
动力 不 能 完全 被 BCOM1 吸收 或 补偿 ””。BDS-3 卫星 的 主体 也 是 长 方 体 ， 因 此 BDS-3 E 
星 定 轨 也 不 适合 继续 用 ECOM1。 


得 比较 


， 因 
验 模 型 。 


3.2.2.3 ECOM2 模型 
时 ， 发 现存 在 伪 谱 线 ， 而 且 GLONASS J 
现 主要 是 由 ECOM1 引起 。 


CODE 分 书 


数据 拟 
此 将 该 先 验 模 型 


在 很 长 时 间 内 ， 对 


日 


由 了 


KE 


[29,30] 


(ECOM1)， 该 模型 在 D 方向 和 YY 方向 仅 保 留 一 个 参 
好 的 定 轨 结 果 ， 这 


D MY 方向 间 的 参数 之 间 


。 该 模型 是 


前 应 用 最 广泛 的 一 个 光 压 模型 


B.cosu+ B, sinu 


im 


29] 


合 得 到 模型 系数 


[入 后 这 种 现象 更 加 


Fi 
bd 


(22) 


期 项 相 比 ，B 方向 的 周期 项 能 显著 地 优化 轨道 模型 的 不 足 。 与 更 早 
新 的 ECOM1， 定 轨 精 度 提 高 有 2~3 倍 。 基 于 ECOM X Block II 和 
4， 确 定 了 6 个 参数 (D、Y、B、2Z1、X1、X3)， 以 及 各 
。 该 模型 精度 相对 于 ROCK 
尺 替 之 前 的 ROCK 模型 ，CODE 从 2005 年 到 
之 后 的 研究 表明 ， 仅 采用 ECOM1 对 GPS 和 GLONASS 进 


高 的 精度 ，2014 年 之 后 的 一 段 时 间 内 CODE 只 是 单独 采 


IT 


HF GPS, GLONASS, BDS-2 等 多 类 卫星 的 光 压 摄 动力 计算 ， 
0D 星 的 定 轨 结果 表明 ， 轨 道 的 激光 检 核 结果 存在 与 8 角 相关 的 系统 


地 球 自 转 参数 和 站 坐标 时 间 序 列 功率 谱 进 


行 分 析 


显 ， 后 来 通 


fF 中 心 的 Arnold 等 人 " (2015) 针对 ECOM1 存在 


步 分 析 发 
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的 问题 ， 进 一 步 对 GPS 和 GLONASS 卫星 进行 研究 ， 结 果 表 明 SRP 仪 产生 沿 D 方向 的 偶 
数 阶 谐 波 的 摄 动 ， 而 对 B 方向 产生 奇数 阶 摄 动 。 该 模型 表达 式 主 要 取决 于 DD M B 方向 周 
期 项 的 数量 ， 通 过 比较 分 析 ， 最 终 从 待定 的 模型 中 选择 了 D 方向 的 两 个 周期 项 和 B 方向 的 
一 个 周期 项 ， 一 共 9 个 参数 ， 建 立 如 下 模型 ， 称 为 BCOM2: 


ap = Do + D2c cos(2Au) + Dos sin(2Au) + D4. cos(4Au) + Das sin(4Au) 
ay = Yo 。 (23) 
ap = Bo + B.cos(Au) + B, sin(Au) 


与 之 前 的 ECOM 系列 模型 不 同 的 是 ， 该 模型 在 D 方向 增加 偶数 项 的 周期 参数 。 此 外 ， 该 模 
型 将 自 变 量 从 久 更 新 到 Au, Au =u- u, 是 轨道 面 上 太阳 与 卫星 之 间 的 轨道 角 ，v。 为 轨道 
面 上 升 交 点 与 太阳 之 间 的 夹 角 。Aw 使 得 周期 参数 的 相位 参考 保持 稳定 ， 与 地 球 公 转 无 关 ， 
可 以 更 方便 地 分 析 估 计 参 数 。 
对 2009—2011 年 的 GPS 和 GLONASS 联合 定 轨 结 果 表 明 ， 与 ECOM1 相 比 ，ECOM2 
可 显著 减少 地 心 坐 标 2 分 量 的 周期 性 误差 ， 而 且 轨 道 天 间 跳 变 、 地 球 自转 参数 的 一 致 性 以 
及 轨道 的 内 符合 精度 也 得 到 提高 。2015 年 1 H4 HAS, CODE 采用 ECOM2 进行 分 析 中 
心 的 产品 计算 。 
虽然 ECOM2 有 其 优越 性 ， 但 是 研究 表明 ECOM2 并 不 是 目前 所 有 卫星 类 型 的 最 优选 
择 。 与 BCOM1 相 比 ， 采 用 ECOM2 的 BDS-2 IGSO 和 MEO 卫星 精密 定 轨 精度 降低 ，。 
对 于 GLONASS 卫星 ， 在 非 地 影 期 时 采用 ECOM1 能 获得 较 高 的 定 轨 精 度 ， 而 在 地 影 区 
时 ， 采 用 ECOM2 的 定 轨 精 度 更 高 ， 这 是 由 于 -X 面 上 的 散热 器 在 地 影 期 产生 扰动 加 速度 ， 
建立 一 个 有 效 的 先 验 模型 可 以 显著 改善 以 上 问题 ”。 对 于 BDS-3 的 SECM 卫星 ， 研 究 发 
现 ， 相 比 于 ECOM1， 基 于 ECOM2 后 激光 检 核 残 差 均 值 更 小 ， 表 明 一 部 分 径 向 误差 可 以 被 
ECOM2 + D 方向 增加 的 周期 项 参数 补偿 ，ECOM2 更 适合 计算 该 类 卫星 的 光 压 摄 动力 ”。 
3.2.3 GSPM 系列 模型 
美国 喷气 动力 实验 室 (Jet Propulsion Laboratory, JPL) 建立 的 GPS 太阳 光 压 模型 
(GPS solar pressure model, GSPM)， 该 模型 用 ( 见 图 1) 的 傅 里 叶 级 数 形式 来 计算 光 压 摄 
动力 ， 也 是 一 种 高 精度 的 光 压 模型 。 
Bar-Sever 在 1997 年 建立 了 GSPM.97 模型 ， 该 模型 在 DY B 坐标 系 中 表示 ， 


D = f(m,r)s(Do + Da cos € + Dsisinet Dez cos 2e + Ds cos 2e) 
Y = Yo + Ya cose + Y,; sine ; (24) 
B= f(m,r)s(Ba cose + By; sine + Beg cos 2e + Bs cos 2¢) 


aH, f(m,r) 是 卫星 质量 和 卫 日 距 x 的 函数 ，s 是 一 个 可 有 效 吸 收 由 于 卫星 本 身材 料及 
环境 变化 等 造成 误差 的 尺度 因子 。GSPM.97 模型 的 待 估 参 数 为 s Do 和 Yoo 

2004 年 ，Bar-Sever 和 Kuang 进一步 开发 了 GSPM.04 模型 ， 采 用 的 坐标 系 为 星 固 
坐标 系 XYZ， 针 对 BlockILA 开发 了 GSPM.ILA.04 模型 ， 针 对 BlockIIR 卫星 开发 了 
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GSPM. II R.04 模型 ， 模 型 参数 为 : 


X = s10~°(AU/r)?/m(X,; sine + Xs sin 2e + Xss sin 3e + Xss sin 5e + Xsr sin 7e) 
| Y = Y + 107-°(AU/r)?/m(Ye1 cos € + Yo cos2e) ; 

Z = 810~-°(AU/r)?/m(Ze1 cos € + Zeg cos 3€ + Zes cos 5é) 

(25) 
其 中 ，AU 是 天 文 单位 ， 模 型 的 待 估 参 数 为 s 和 Yh。 新 的 模型 与 之 前 模型 相 比 ，Block IIA 
与 BlockIIR 卫星 的 轨道 拟 合 精度 分 别提 高 了 24% 和 80%， 轨 道 预报 精度 分 别提 高 了 32% 
和 58 名。 在 此 基础 上 ， 分 别 采 用 1997 年 1 月 到 2010 年 5 月 和 1992 年 到 2013 年 的 数据 ， 
根据 卫星 类 型 进一步 开发 了 新 版 本 GSPM10 和 GSPM13””。 

JPL 采用 GIPSY-OASIS 软件 的 研究 结果 表明 ， 基 于 GSPM 的 方法 对 于 轨道 重 登 弧 段 
精度 、 激 光 检 核 残 差 、 站 坐标 重复 性 等 方面 可 以 得 到 更 精确 的 解 ， 而 ECOM1 SF RBM 
段 钟 差 和 日 长 变化 能 得 到 更 高 的 精度 。 在 2008 年 时 ，8 个 IGS 分 析 中 心中 的 5/6 个 采用 的 
为 ECOM1， 其 中 有 1/2 个 没有 采用 先 验 模型 ”。 对 于 GPS Block IF 和 IIR PÆ, IGS 第 
三 次 重 处 理 推 荐 了 ECOM1， 除 此 之 外 的 其 他 GPS 卫星 ， 推 荐 采用 ECOM2 或 GSPM。 对 
F BDS-3, 与 ECOM2 相 比 ，Box-Wing 模型 作为 先 验 模型 的 ECOMI1 能 获得 更 高 的 轨道 
精度 ""。 对 于 经 验 模 型 ， 不 顾及 二 阶 项 的 七 参数 ECOM? 模型 能 获得 较 好 的 BDS-3 卫星 轨 
道 精 度 ， 其 中 重 闪 轨道 精度 和 SLR 检 核 精 度 分 别 约 为 6 cm 和 4 cm". 

与 分 析 型 模型 相 比 ， 经 验 型 光 压 模型 参数 与 轨道 其 他 参数 同时 估计 ， 不 再 依赖 于 卫 
的 光 压 参数 ， 而 且 表 达 式 更 简洁 ， 精 度 更 高 。 经 验光 压 模型 是 影响 轨道 精度 的 主要 因素 ”。 
但 是 经 验 型 模型 也 有 一 定 的 缺点 : (1) 该 类 模型 不 是 对 卫星 受到 的 光 压 摄 动力 进行 受 力 分 析 
后 建立 的 ，(2) 可 能 会 与 其 他 的 非 保守 力 耦 合 ， 在 太阳 高 度 角 较 大 时 可 能 会 与 卫星 的 相位 中 
心 偏差 等 参数 相关 ; (3) 同一 个 经 验 模型 并 不 具有 普 适 性 ， 对 不 同类 型 的 卫星 需要 建立 并 采 
用 不 同 的 经 验 模型 。 
3.3 ”半分 析 型 光 压 模型 
3.3.1 可 校正 Box-Wing 模型 

由 于 Box-Wing 模型 需要 卫星 的 准确 尺寸 以 及 光 压 参数 等 ， 在 没有 这 些 参数 的 准确 信 
息 时 ， 难 以 应 用 Box-Wing 模型 ”。Rodriguez-Solano 等 人 (2012) 开发 了 可 以 直接 估计 卫 
星光 压 参数 的 可 校正 Box-wing 模型 (adjustable Box-Wing, ABW)， 该 模型 属于 半分 析 模 
测 户 

民 据 式 (17) 和 (18) 可 以 计算 太阳 能 帆 板 和 卫星 本 体 受到 的 光 压 力 。 基 于 名 义 姿 态 模型 
和 角度 =， 可 以 得 到 加 速度 对 各 光 压 参数 的 偏 导 数 ， 进 而 对 以 上 参数 进行 估计 。 对 于 太 虹 
能 帆 板 ， 其 受到 的 光 压力 在 日 国 华 标 系 中 表示 ，ep 与 ex 的 方向 一 致 ，cos9 = 1， 根 据 式 
(17)， 则 其 光 压 参数 只 有 一 个 为 1+p 十 35; 对 于 卫星 本 体 的 受 照 面 ， 其 受到 的 光 压力 在 星 
固 坐 标 系 中 表示 ， 根 据 式 (18), +X. +2 -Z 三 个 受 照 面 的 两 个 光 压 参数 为 a 十 6 和 p, 
£6 个 参数 。ABW 模型 估计 的 参数 共 9 个 ， 另 外 2 个 为 轴 Y 偏差 参数 (太阳 能 帆 板 的 实际 
对 准 角 度 与 设计 存在 微小 偏差 ) 以 及 卫星 帆 板 的 旋转 滞后 角 参 数 。 
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J 计算 的 轨道 精度 与 BCOMI1 相当 ， 但 
目前 该 模型 并 
模型 较 复杂 ， 各 参数 之 间 的 相关 性 较 强 ， 


是 基于 两 个 模型 的 轨道 之 间 有 系统 差 ， 
未 被 应 用 于 各 IGS 分 析 中 心 日 常 的 定 
模型 精度 依赖 于 先 验 卫星 参数 


和 合理 约束 。 
3.3.2” 先 验 模型 的 构建 

HF EA 
际 的 物理 意义 ， 常 用 的 策略 为 一 个 
近年 来 关于 光 压 模型 的 研究 很 


因此 ， 该 模型 主要 对 了 


F 光 压 参 数 未 知 


Ee: 


EY 其 先 验 光 压 模型 。 
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E 验 模型 并 以 经 验 模型 参数 吸收 未 模型 化 的 误差 ， 因 
要 的 一 个 方向 就 是 先 验 模型 的 构建 。 
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目前 比较 常用 的 2 


模型 为 Box-Wing 模型 ， 与 单独 的 ECOM1 和 ECOM2 模型 相 比 ， 该 模型 应 用 于 Galileo, 


QZSS、BDS-3e、BDS-3 等 非 立 方 体 结构 ] 
几何 参数 和 光学 特性 参数 "”"。 除 此 之 外 ， 有 些 学 者 也 提出 了 多 种 
对 Galileo-IOV 卫星 在 采用 ECOMI1 i 


差 ， 该 误差 与 卫星 的 太阳 高 度 角 相关 ， 


压 摄 动 力 未 被 光 压 模型 有 效 吸收 造成 ”。 
Raumfahrt, DLR) 的 Montenbruck 等 人 ”基于 


(DLR 模型 )， 该 模型 从 最 原 
部 分 : 立方 体 、 伸 长 、 非 对 称 结构 ， 
动 加 速度 (Y 方向 在 名 义 姿 态 下 不 受 光 
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来 吸收 )， 
的 先 验 值 。DLR , 模型 做 先 验 模型 的 ECOMI1 可 消除 上 述 周 期 特 
主体 为 非 正方 体 的 卫星 的 太阳 光 压 摄 动力 。 


始 的 Box-Wing 模型 
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得到 了 卫星 本 体 
压 摄 动力 影响 
玫 抽 象 为 少量 与 每 部 分 结构 的 a 二 6 和 p 相关 的 参数 的 表达 式 ， 
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4 类 似 的 方法 ， 选 择 了 6 Set RA 
验光 压 模 型 参数 估 值 。 


进行 精密 定 轨 获 得 精密 光 压 摄 动 加 速度 ， 以 精密 光 压 摄 动 加 


上 
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发 ， 将 卫星 主体 受到 的 加 速度 分 为 三 个 
在 D、B (或 者 X、2) 方向 的 光 压 摄 
， 如 果 受 到 光 压 力 ， 可 以 通过 ECOM1 
并 估计 了 该 参数 
日 于 计 


误差 ， 表 明 该 模型 适 上 有 


型 参数 ， 采 用 长 期 的 
基体 实现 为 :， 先 用 ABW 模型 
速度 为 观测 值 平 差 获 取 6 个 先 


验光 压 模 型 参数 ， 每 天 估计 出 一 组 参数 ， 
参数 ， 根 据 计 算 的 参数 ， 得 到 光 


压 关 于 s 的 函数 对 
精密 定 轨 精 度 并 减弱 钟 差 拟 合 残 差 随 轨道 的 系统 必 
Wang 等 人 "考虑 综合 ECOM1 和 ABW 模型 (顾及 零 偏 姿态 、 通 


最 终 获 得 


一 组 最 终 的 参数 ， 作 为 


固定 的 ABW 模型 


长 达 式 。 结 果 表 明 ， 该 先 验 模 型 能 显著 提升 
FE 变化 误差 。 


的 优点 ， 构 建 了 BDS-2 GEO 卫星 的 光 


压 模 型 。 
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后 者 应 用 于 零 偏 模式 下 ); 对 不 同方 向 
项 式 ， 为 了 减少 参数 ， 分 别 对 和 轨道 
分 析 (选择 影响 较 大 的 角度 和 项 ， 在 D 
虑 正弦 项 ); 基于 ABW 计算 的 SRP 对 


为 自 变 量 的 二 阶 多 项 式 。 厦 


dade 


角 以 及 


采用 ABW 模型 的 单 天 解 计 
帆 板 与 太阳 光 并 不 垂直 ， 因 此 DYB } 
示 零 偏 下 的 太阳 光 压 ， 因 此 将 其 分 解 到 DYB M DYB 方向 (前 者 应 
的 光 压 进行 谱 


言 天 线 及 谈 挡 影响 ) 
SRP; 由 于 
完全 适用 于 


/| F% 
航模 式 下 ， 


E 架 并 不 
IE SU 
分 析 ， 获 得 其 周期 项 ， 建 立 三 角 函 数 多 
Ea AASB eA AG EE A TLR 


BRD AALM SRR, FEY MB 方向 仅 考 


参数 进行 拟 合 ， 分 别 表示 为 线性 模型 、 大 PB 
究 结果 表明 ， 该 模型 能 够 显著 消除 


高 度 角 GB 
日 太阳 光 压 导致 的 与 太阳 方位 


角 、 轨 道 角 有 关 的 轨 


da 


误差 ， 降 低 轨 道 径 向 的 系统 偏差 。 
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图 4 为 自 导 航 卫 星 发 射 以 来 光 压 模型 研究 的 发 展 历程 。 基 于 以 上 分 析 可 以 发 现 ， 目 前 


常用 的 导航 卫星 定 轨 的 光 压 模型 较 多 ， 各 种 模型 的 适用 性 主要 与 卫星 形状 和 设计 有 关 。 表 2 
列 出 目前 导航 卫星 常用 光 压 模型 的 一 些 信息 ， 其 中 适用 卫星 是 已 有 文献 的 分 析 结 果 ， 对 于 一 
些 新 卫星 适用 的 光 压 模型 还 需要 进一步 研究 论证 。 


es 
ROCK T10/T20 
1999 
ogre ECOM1 
1985 ROCK42 \ _ [1996 >` 
ROCK4 Colombo \ ROCK T30 
| | = = = a = 
1980 2020 
1994 2004 2010 2013 
1997 
ECOM cspMo7 GSPM.04 GSPM.10 GSPM.13 
4 ”导航 卫星 光 压 模型 发 展 时 间 轴 
表 2 目前 导航 卫星 的 推荐 光 压 模型 
ECOM 系列 模型 ECOM1+Box 
费 型 名 称 GSPM 模型 ABW 
aa ECOMI ECOM2 Re Wing 
坐标 系 DYB DYB XYZ XYZ DYB 
模型 类 别 经 验 型 经 验 型 经 验 型 半分 析 型 半分 析 型 
参数 个 数 5 9 2 9 5 
GPS-IIR-M IIA 
GPS-IIR-M IIIA 
GLONASS GPS QZSS 
Qe GPS IIA IIR IIF GPS-IIA IIR QZSS 
适用 卫星 QZSS Galileo 
BDS-2 IIR-M IIIA Galileo 
Galileo BDS 
BDS 
BDS-3 
4 4 论 


高 精度 的 卫星 轨道 是 实现 导航 卫星 高 精度 应 用 的 前 提 ， 只 有 精细 地 考虑 卫星 受到 的 各 
' 摄 动力 才能 精确 确定 卫星 轨道 ，SRP 作为 卫星 所 受 非 保守 力 的 主要 来 源 需 要 精密 建 模 。 
姿态 模型 是 建立 光 压 模型 的 基础 ， 姿 态 模型 准确 与 否 会 直接 影响 到 光 压 模型 精度 。 本 文 介 绍 
了 目前 主要 导航 卫星 的 姿态 模型 和 常用 的 光 压 模型 ， 详 细 分 析 了 有 具体 模型 参数 以 及 各 自 的 
应 用 ， 总 结 比较 了 不 同 光 压 模 型 的 优 缺 点 。 
目前 各 分 析 中 心 常 用 的 光 压 模型 为 ECOM 系列 模型 或 者 GSPM 系列 模型 ， 并 采用 相 
应 的 先 验 模型 ， 这 两 个 系列 模型 基本 能 满足 目前 各 类 导航 卫星 的 精密 定 轨 需求 。 但 是 目前 
关于 导航 卫星 姿态 和 光 压 模型 还 需要 进一步 研究 ， 部 分 BDS-3 卫星 的 实际 姿态 与 目前 公 
布 的 姿态 模型 并 不 完全 一 致 ， 轨 道 机 动 时 的 姿态 变化 需要 进一步 研究 ， 由 于 缺乏 卫星 的 详 


细 结 构 参 数 ， 对 于 卫星 在 其 运行 不 同时 期 尤其 是 地 影 期 的 光 压 建 模 精 度 较 低 ， 由 于 不 同 卫 
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星 的 结构 设计 有 差别 ， 建 立 分 析 型 模型 时 需要 区 分 考虑 ， 这 需要 对 不 同 卫 星 独 立 分 析 ， 例 
W BDS GEO 卫星 的 大 天 线 会 对 卫星 的 光 压 有 较 大 影响 ， 目 前 缺少 与 最 新 发 射 的 GPS IR 
GLONASS-K 和 BDS-3 IGSO 卫星 的 光 压 模型 相关 的 文献 ， 需 要 进一步 研究 。 
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Research Progress of Attitude and Solar Radiation Pressure 


Model for Navigation Satellites 


HUANG Chao!?, SONG Shu-li!, CHEN Qin-ming!, ZHOU Wei-li!? 


(1. Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai 200080, China; 2 
University of Chinese Academy of Sciences, Beijing 100049, China) 


Abstract: The shape, material, orbit design, and attitude control model are different for 
satellites of different navigation systems. In order to obtain high-precision satellite orbits, 
different types of satellites need to be processed differently in terms of satellite attitude 
control mode and solar pressure model. In order to establish the best model to obtain higher 
orbit determination accuracy, there have been a lot of works focused on the establishment 
and refinement of the solar pressure model and attitude control mode of navigation satellites 
because solar pressure is the largest nonconservative force. Recent research shows that 
the new Galileo and Beidou satellites have adopted different attitude control modes and 
solar radiation pressure models. This paper summarizes the general attitude models of 
navigation satellites, especially the attitude model of Beidou satellite, which is officially 
published recently. The research progress of the solar radiation pressure model is introduced 
and the advantages and disadvantages of each model and their adaptability are analyzed. 
Through the analysis of the current attitude and solar radiation pressure models, the needed 
improvement in higher precision orbit determination is proposed. This paper provides a 
reference for the selection and establishment of solar radiation pressure and attitude models 
in the precision orbit determination of navigation satellites, especially for newly launched 


satellites. 


Key words: GNSS satellite; solar radiation pressure; satellite attitude; precise orbit deter- 


mination 


